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In this talk, I start with a brief introduction to Non-Relativistic QCD (NRQCD) and 
^ . its apphcations to quarkonium physics. This theory has provided a consistent frame- 

work for the physics of quarkonia, in particular, the colour-octet Fock components 
predicted by NRQCD have important implications for the phenomenology of charmo- 
nium production in experiments. The applications of NRQCD to J/ip production at 
Tevatron and the tests of the theory in other experiments is discussed. In particular, 
the apparent disagreement of NRQCD with results from HERA on inelastic photopro- 
duction of J/ip is discussed and it is shown that the results are rather susceptible to 
intrinsic transverse momentum smearing. The photoproduction data, therefore, do not 
provide a good test of NRQCD. It is argued that NRQCD may be tested stringently 
by looking for the production of other charmonium resonances at the Tevatron, be- 
cause the production rates for these resonances can be predicted within the NRQCD 
framework. 



ABSTRACT 



*sridhar@theory.tifr.res.in 



Over the last few years, there has been a considerable advance in the understanding 
of quarkonium physics due to the development of the non-relativistic effective field 
theory of QCD, called non-relativistic QCD (NRQCD) [l[]. The Lagrangian for this 
effective theory is obtained from the full QCD Lagrangian by neglecting all states with 
momenta larger than a cutoff of the order of the heavy quark mass, m, and accounting 
for this exclusion by introducing new interactions in the effective Lagrangian, which 
are local since the excluded states are relativistic. Beyond the leading order in l/m 
the effective theory is non-renormalisable. The scale m is an ultraviolet cut-off for the 
physics of the bound state; however the latter is more intimately tied to the scales 
mv and mv^, where v is the relative velocity of the quarks in the bound state. The 
physical quarkonium state admits of a Fock expansion in v, and it turns out that the 
QQ states appear in either colour- singlet or colour-octet configurations in this series. 
Of course the physical state must be a colour- singlet, so that a colour-octet QQ state 
is connected to the physical state by the emission of one or more soft gluons. In 
spite of the non-perturbative nature of the soft gluon emissions, the effective theory 
still gives useful information about the intermediate octet states. This is because the 
dominant transitions that occur from colour-octet to physical colour-singlet states are 
via El or Ml transitions with higher multipoles being suppressed by powers of v. It 
then becomes possible to use the usual selection rules for these radiative transitions 
to keep account of the quantum numbers of the octet states, so that the production 
of a QQ pair in a octet state can be calculated and its transition to a physical singlet 
state can be specified by a non-perturbative matrix element. The cross-section for the 
production of a meson H then takes on the following factorised form: 

^(H) = Y. -^^.iO^C'-'Lj)), (1) 

n={a,S,L,J} '"" 

where the F„'s are the short-distance coefficients and the On are local 4-fermion oper- 
ators, of naive dimension dn-, describing the long-distance physics. The short- distance 
coefficients are associated with the production of a QQ pair with the colour and angu- 
lar momentum quantum numbers indexed by n. These involve momenta of the order 
of m or larger and can be calculated in a perturbation expansion in the QCD cou- 
pling as(m). The QQ pair so produced has a separation of the order of l/m which 
is pointlike on the scale of the quarkonium wavefunction, which is of order l/{mv). 
The non-perturbative long-distance factor {O^) is proportional to the probability for 
a pointlike QQ pair in the state n to form a bound state H . 

The existence of the colour-octet components of the quarkonium wave function is 
the new feature of the NRQCD approach. Before the development of NRQCD, the 
production and decay of quarkonia were treated within the framework of the colour- 
singlet model 0, 1^. In this model, it is assumed that the QQ pair is formed in the 
short- distance process in a colour-singlet state. The corrections from terms higher 
order in v were neglected. While this model gave a reasonable description of low- 
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energy J /ip data, it was known that it was incomplete because of an inconsistency in 
the treatment of the P-state quarkonia. This was due to a non-factorising infra-red 
divergence, noted first in the apphcation of the colour- singlet model to Xc decays |Q, 
and the proper resolution of this problem was obtained only by including the colour- 
octet components in the treatment of the P-states 0. The colour-octet components, 
however, had a more dramatic impact [^ on the phenomenology of P-state charmonium 
production at large pt at the Tevatron pp collider 0] where the colour-singlet model 
was seen to fail miserably. While the inclusion of the colour-octet components for the 
P-states was necessary from the requirement of theoretical consistency, there was no 
such problem with the S states because the corresponding amplitude was finite and the 
colour-octet components were suppressed compared to the colour-singlet component 
by 0{y^). But the data on direct J/ip and ip' production at the Tevatron [0 seem to 
indicate an important contribution from the colour-octet components for the S-states 
as well §. 

While it is clear that the correct description of the Tevatron large-pr data re- 
quires that the colour-octet components of the quarkonium wave function have to be 
taken into account, the major problem is that the corresponding long-distance matrix 
elements are a priori unknown and can be obtained only by fitting to the Tevatron 
data 1^, |10[- The direct J/ip production cross section in the NRQCD approach re- 



ceives contributions from the colour-singlet ^5*1 channel and the colour-octet '^Pj\ 
^Sq and ^5*1 channels. The non-perturbative parameter for the colour-singlet chan- 
nel is known from J /if) leptonic decay. Given this input, the three non-perturbative 
parameters {Oip^Pf)), {0{^Sf^)), {0{^Sf^)) (which we call matrix elements Mi, Ms 
and M3 respectively) are extracted from a fit to the CDF data. It turns out that for 
Pt > 4: GeV, the pt dependence of the short- distance coefficients corresponding to the 
^Pj and the ^Sq channels are identical. The ^S[ channel on the other hand has 
a different pt distribution, because fragmentation-type contributions are present only 
only for this channel. Consequently, the shape of the experimental pt distributions 
can be used to determine M3 spearately, but only a linear combination of Mi and M2 
(i.e. Ml/ml + M2/3) can be fitted. 

Clearly it is important to have other tests of NRQCD, and much effort has been 
made recently to understand the implications of these colour-octet channels for J/ip 
production in other experiments. We discuss some of these below. 



1. The prediction |[TT| , [T2|, [T^ for prompt J/ip production at LEP in the colour- 
singlet model is of the order of 3 x 10~^, which is almost an order of magnitude 
below the experimental number for the branching fraction obtained from LEP 



T^. Recently, the colour-octet contributions to J/^p production in this chan- 



nel have been studied |T5|, |I^ and it is found that the inclusion of the colour- 
octet contributions in the fragmentation functions results in a predictions for the 
branching ratio which is 1.4 x 10~^ which is compatible with the measured values 



of the branching fraction from LEP IT^. A more accurate analysis, resumming 



large logarithms in Ejj^/Mz ignored in Ref. |0 has been recently performed 



2. The production of J/ip in low energy e^e machines can also provide a stringent 
test of the colour-octet mechanism |TB[. In this case, the colour-octet contri- 



butions dominate near the upper endpoint of the J/ip energy spectrum, and 
the signature for the colour-octet process is a dramatic change in the angular 
distribution of the J/^ near the endpoint. 

One striking prediction of the colour-octet fragmentation process both for pp 
colliders and for J/ip production at the Z-peak, is that the J /if) coming from the 
process g — > J/ipX is produced in a transversely polarised state JTOf. For the 



colour-octet cc production, this is predicted to be a 100% transverse polarisation, 
and heavy-quark spin symmetry will then ensure that non-perturbative effects 
which convert the cc to a J/ip will change this polarisation only very mildly. This 
spin-alignment can, therefore, be used as a test of colour-octet fragmentation. 



4. The colour-octet components are found [Q to dominate the production processes 
in fixed-target pp and vrp experiments. Using the colour-octet matrix elements 
extracted from elastic photoproduction data it is possible to get a very good de- 
scription of the -N/i-dependence and also the xp and rapidity distributions. More 
recently, NLO corrections to the fixed-target cross-sections have been calculated 
121. 



5. The associated production of a J/ip + 7 is also a crucial test of the colour-octet 
components |^ and also of the fragmentation picture [^. Similar tests can be 
concieved of with double J /ip production at the Tevatron ^A . 



6. J/ip and ip' production in pp collisions at centre-of-mass energies of 14 TeV at 
the LHC also provides a crucial test of colour-octet fragmentation ||2^. Recently, 
J/ if) + 7 production at the LHC has also been studied [^ at the LHC. 



One important cross-check is the inelastic photoproduction of J/ip at the HERA 
ep collider |^. The inelasticity of the events is ensured by choosing z = Pp-pj/^/Pp-p-y 



to be sufficiently smaller than one and, in addition, using px > I GeV. The surprising 
feature of the comparisons [^, ^ of the NRQCD results with the data from HERA is 
that the colour-singlet model prediction is in agreement with the data while including 
the colour-octet component leads to violent disagreement with the data at large z. 
While the colour-singlet cross section dominates in most of the low- 2; region, the colour- 
octet contribution increases steeply in the large- 2; (0.8 < z < 0.9) region and this rise is 
not seen in the data. In these comparisons, the values of the non-perturbative matrix 
elements are taken to be those determined from a fit to the Tevatron large-pr data. 



Naively, one would think that this points to a failure of NRQCD. But this conclusion 
is premature. The reason is that while at the Tevatron the measured pt of the J/ip 
is greater than about 5 GeV, at HERA the pr can be as small as 0{1) GeV. At 
such small values of pt (and also for z very close to unity), there could be significant 
perturbative and non-perturbative soft physics effects. One way to explore the effect of 
such contributions is to include transverse momentum smearing of the partons inside 
the proton, and study the effects of the parton transverse momentum, k^, on the 



J/ijj distributions both at the Tevatron and at HERA. It has been demonstrated |30] 
that the z distribution measured at HERA is particularly sensitive to the effects of kx 
smearing, and that inelastic photoproduction at HERA, with the present kinematic 
cuts, is not a clean test of NRQCD [} 

In Fig. 1, the results of the fits to the Tevatron data are shown, for three different 
values of average kx, {kr), viz. (kr) = 0, 0.7, 1.0 GeV. It is observed that the effect of 
the kx smearing on the parameters extracted from the data is very modest. Fig. 1 shows 
that the fits to the data when kx smearing is included are very good and comparable in 
quality to the case (kx) = 0. Taking these fitted values of the parameters, inelastic J/ip 
photoproduction at HERA is considered, for the same choice of parton distributions, 
scales etc. as used in the Tevatron fits. The z distribution, for ^Js^p = 100 GeV and 
Pt > 1 GeV, is compared with the data from HERA in Fig. 2. Again the theoretical 
curves in Fig. 2 are for (kr) = 0, 0.7, 1.0 GeV. 

In the absence of smearing, (fcy) = 0, we see that the colour-octet component makes 
a large contribution at z close to 1 which is not supported by the data. However the 
introduction of kj- makes a substantial change to the octet contribution. Whereas the 
effect of fcr-smearing is very small for large-pr production at the Tevatron, these effects 
are found to be very important for J/ip production at HERA. In particular, smearing 
signficantly reduces the size of the cross section and the z distribution also becomes 
flatter, in better agreement with the HERA data. It is safe to conclude that while a 
direct comparison of the NRQCD predictions with the ^-dependence of the inelastic 
photoproduction cross section for J/ip at HERA show a marked disagreement between 
the two, we argue that such a comparison is misleading. The inelastic photoproduction 
process does not provide a clean test of NRQCD because of the very low px-cut (~ 
1 GeV) used in the HERA experiments making the data very susceptible to effects like 
kx smearing. 

Better tests of NRQCD may be obtained by studying other observables at the 
Tevatron itself. The study of the polarisation of the produced J/ip mentioned earlier 
is one example; in the following, we will discuss the production of other charmonium 
resonances [p3| , whose cross-section can be predicted in NRQCD. One important 
feature of the NRQCD Lagrangian is that it shows an approximate heavy-quark sym- 
metry, which is valid to 0{v^) ~ 0.3. The implication of this symmetry is that the 

^ Other effects such as soft-gluon resummation |gl[ and the breakdown of NRQCD factorisation 
near z = 1 [p2| have been discussed in the context of this discrepancy. 
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Figure 1: The CDF data for Bda/dpT (in nb/GeV) for J/-?/; production at 1.8 TeV 
with —0.6 < T] < 0.6, compared to the model predictions with (kr) = 0, 0.7, 1.0 GeV, 
respectively. 



nonperturbative parameters have a weak dependence on the magnetic quantum num- 
ber. Using this symmetry some non-perturbative matrix elements can be expressed in 
terms of others already determined from the Tevatron data. In particular, the ^Pi ma- 
trix elements can be inferred from the Tevatron data on x production and, therefore, 
the production of the ^Pi charmonium state, he, can be predicted in NRQCD |^ . 
The production of the he is interesting in its own right : charmonium spectroscopy 
predicts this state to exist at the centre-of-gravity of the xd^Pj) states. While 



the E760 collaboration at the Fermilab has reported ||3^ the first observation of this 
resonance its existence needs further confirmation. 

The cross-section for he production at the Tevatron energy (i/i = 1.8 TeV) has 
been presented in Ref. [^ For 20 pb~^ total luminosity, for pt integrated between 5 
and 20 GeV we expect of the order of 650 events in the J/ip + ■n channel. Of these, the 
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Figure 2: The HERA data [^ for da/dz (in nb) for J /ip production at J's^ = 100 GeV 
with pt > I GeV, compared with model predictions for three choices of the intrinsic 
transverse momentum distribution, namely (/ct) = 0, 0.7, 1.0 GeV. 



contribution from the colour-singlet channel is a little more than 40, while the octet 
channel gives more than 600 events. The colour-octet dominance is more pronounced at 
large-pr- Recent results on J/ip production from CDF are based on a total luminosity 
of 110 pb~^. For this sample, more than 3000 events can be expected to come from the 
decay of the he into a J/ip and a tt. With this large event rate, the he should certainly 
be observable if the pi° coming from its decay can be reconstructed efficiently. 

A similar prediction for the absolute production rate can be made for r/c production 
where the two-photon decay mode of the rjc has been considered. 
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symmetry allows the 7]^ cross-section to be determined in terms of the non-perturbative 
parameters Mi and M2 obtained from J/ijj data. But as explained before, the J/ip 
data do not allow for a separate determination of these parameters, but only a linear 
combination of these parameters. We can saturate the linear combination with either 



Ml or M2 and we obtain the rjc event rate in both these cases. For the integrated event 
rate, with a pT-cut of 5 GeV and assuming an integrated luminosity of 110 pb~^, we 
find that the number of r^c — > 77 hes between 425 and 7700, depending on whether Mi 
or M2 saturates the hnear combination. The sensitivity of the event rate to Mi and 
M2 shows that the experimental measurement of the etac cross-section will allow for 
an accurate determiantion of these non-perturbative parameters. We reiterate that the 
rates for he and rjc are predictions of NRQCD, and it is not possible to have similar 
predictions in alternative approaches to quarkonium production like colour-evaporation 

In conclusion, NRQCD provides a predictive theoretical framework for quarko- 
nium physics. In particular, the anomalies in the J/ip production at the Tevatron are 
properly understood using NRQCD. Several other tests of the theory, proposed in the 
literature, have been discussed. In particular the inelastic photoproduction of J/ip at 
HERA is discussed and it shown that the apparent disagrement of the experimental 
results with the predictions of NRQCD is misleading. Because of the low values of pt 
in the photoproduction case, we find that the effect of fc-r-smearing is important and 
that, indeed, for (kx) ~ 0.7 GeV, the discrepancy between theory and experiment is 
no longer observed. On the other hand, the inclusion of kj- smearing has a very mod- 
est effect on the large-pr J/^^ data from the Tevatron. Better tests of NRQCD may 
be obtained by studying other observables at the Tevatron itself, such as the study 
of the polarisation of the produced J/ip |T^ or the production of other charmonium 
resonances |^, ^ whose cross-section can be predicted in NRQCD. 
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